Bosonization, Pairing, and Superconductivity of the Fermionic Tonks-Girardeau Gas 
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We determine some exact static and time-dependent properties of the fermionic Tonks-Girardeau 
(FTG) gas, a spin-aligned one-dimensional Fermi gas with infinitely strongly attractive zero-range 
odd-wave interactions. We show that the two-particle reduced density matrix exhibits maximal 
superconductive off-diagonal long-range order, and on a ring an FTG gas with an even number of 
atoms has a highly degenerate ground state with quantization of Coriolis rotational flux and high 
sensitivity to rotation and to external fields and accelerations. For a gas initially under harmonic 
confinement we show that during an expansion the momentum distribution undergoes a "dynamical 
bosonization" , approaching that of an ideal Bose gas without violating the Pauli exclusion principle. 
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If an ultracold atomic vapor is confined in a de Broglie 
wave guide with transverse trapping so tight and tem- 
perature so low that the transverse vibrational excita- 
tion quantum Hlo is larger than available longitudinal 
zero point and thermal energies, the effective dynam- 
ics becomes one-dimensional (ID) 0, 0j, a regime cur- 
rently under intense experimental investigation [3, 
Confinement-induced ID Feshbach resonances (CIRs) 
reachable by tuning the ID coupling constant via 3D Fes- 
hbach scattering resonances occur for both Bose gases flj 
and spin-aligned Fermi gases |^. Near a CIR the ID in- 
teraction is very strong, leading to strong short-range cor- 
relations, breakdown of effective-field theories, and emer- 
gence of highly-correlated iV-body ground states. In the 
bosonic case with very strong repulsion (ID hard-core 
Bose gas with coupling constant +oo, the Tonks- 

Girardeau (TG) gas), the exact iV-body ground state 
was determined some 45 years ago by a Fermi-Bose (FB) 
mapping to an ideal Fermi gas [6j, leading to "fermion- 
ization" of many properties of this Bose system, as re- 
cently confirmed experimentally Q . The "fermionic TG" 
(FTG) gas 0, a spin-aligned Fermi gas with very strong 
attractive ID odd-wave interactions, can be realized by 
3D Feshbach resonance mediated tuning to the attractive 
side of the CIR with ID coupling constant g[jj — oo. It 



has beenpointed out [a, lD| that the generalized FB map- 
ping lajldlal can be exploited in the opposite direction to 
map this system to the trapped ideal Bose gas, leading 
to determination of the exact iV-body ground state and 
"bosonization" of many properties of this Fermi system. 
We recently examined the equilibrium one-body density 
matrix and exact dynamics following sudden turnoff of 
the interactions by detuning from the CIR 0. Here we 
determine some other exact properties of the untrapped, 
ring-trapped, and harmonically trapped fermionic TG 
gas, the most striking of which are pairing, superconduc- 
tive off-diagonal long-range order (ODLRO) of the two- 
body density matrix, a highly degenerate ground state 



of an even number of atoms on a ring with quantiza- 
tion of Coriolis rotational flux and high sensitivity to 
rotation and to external fields and accelerations, and a 
"dynamical bosonization" of the momentum distribution 
following sudden relaxation of the trap frequency. 

Untrapped FTG gas: The Hamiltonian is H = 

'Ef=i [-|^^J + J2i<j<e<N vfnti^j - ^e) where <t is 
the two-body interaction. Since the spatial wave func- 
tion is antisymmetric due to spin polarization, there 
is no s-wave interaction, but it has been shown |^ Q 
that a strong, attractive, short-range odd-wave interac- 
tion (ID analog of 3D p-wave interactions) occurs near 
the CIR. This can be modeled by a narrow and deep 
square well of depth Vq and width 2xn. The contact con- 
dition at the edges of the well is Q ■tppixji ~ xq) ^ 
-'0f(2^j£ = -a;o) = -afj~,ijj'p{xji = ±Xo) where af^ is 
the ID scattering length and the prime denotes differen- 
tiation. Consider first the relative wave function ippix) 
in the case N = 2. The FTG limit is af^ — oo, a 
zero-energy scattering resonance. The exterior solution 
is iPf{x) = sgn(a;) = ±1 (-1-1 for a; > and —1 for x < 0) 
and the interior solution fitting smoothly onto this is 
sin(Kx) with k = mVo/h? = t:/2xo. In the zero-range 
limit xq 0+ the well area 2x0 Vb = {ttH)'^ /2mxo oo, 
stronger than a negative delta function. In this limit the 
wave function is discontinuous at contact xq = 0±, al- 
lowing an infinitely strong zero-range interaction in spite 
of the antisymmetry of ipF 0- This generalizes immedi- 
ately to arbitrary A^: the exact FTG gas ground state is 



N 



IppfyXl 



,Xn) = A{xi,- ■ ■ ,XN)W_<i)o{Xj) (1) 



with A{xi, ■ ■ ■ ,xn) = ni<j<£<Ar sgn(x£ - Xj) the "unit 
antisymmetric function" employed in the original discov- 
ery of fermionization and 0o = 1/ the ideal Bose 
gas ground orbital, L being the periodicity length. Its 
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energy is zero [lOj and it satisfies periodic boundary con- 
ditions for odd A'^ and antiperiodic boundary conditions 
for even TV [ill- 

The exact single-particle density matrix pi{x,x') — 

N r jlf'jx, X-2, - ■ ■ ,XN)tp*F{x' ,X2, ■ ■ ■ ,XN)dX2 - ■ ■ dXN 

isHHa Pi{x.x') = N(l)o{x)(t)*Q{x')[F{x,x')]^-^ with 
F(x,x') = /_[/2Sgn(j; - y)sgn(a;' - y)\My)?dy = 
\ — 2\x — x'\/ L. In the thermodynamic limit N ^ oo, 
L — !■ oo, N/L — n this gives an exponential decay 
[l^ : pi{x,x') = ne~^"l^^^ Its Fourier transform Uk, 
normalized to Uk = N (allowed momenta viii / L 
with V = Q, ±1, ±2, • • • ), is the momentum distribution 
fmiction n/j = [1 + {k/2nY]^^ . It satisfies the exclusion 
principle limitation rtfe < 1, but nevertheless, for n — > 
the continuous momentum density nik) = (L/27r)nfc 
reduces to N times a representation of the Dirac delta 
function, simulating the ideal Bose gas distribution: 
n{k);^N5{k) [13. 

The two-particle density matrix p2{xi,X2\x'i^x'2) — 

N{N -1)/ iPf{xi, ■■,XN)-4'*p{x''i^,x'2,Xs„ ..,XN)dxs,...dxN 

also has a simple closed form: 

P2{xi,X2]x[,x'2) = N{N - l)sgn(a;i - X2)(j)o{xi)(j)o{x2) 
X sgn(i;-4)0*(x;)0S(4)[G(a;i,a;2;a;;,4)]~"2 (2) 

where [G(a:i, 0:2; xi, X2)]^"^ = [/_^^^2 '^gn(a;i - 

a::)sgn(a;2 — x)sgn{x'i — a;)sgn(a;2 — x)\(j)f)\^ {x)dx\^ "'^ = 
g2n(yi-y2+y3-yi) ^^i^ thermodynamic limit and 

yi < 2/2 < 2/3 < 2/4 are the arguments {xi,X2]x'i,x'2) 
in ascending order. p2 is of order in the following 
cases: (a) \xi — x'i\ < 0(l/n), \x2 — x'2\ < 0{l/n); 
(b) \xi - 41 < 0(l/n), \X2 - x[\ < 0(l/n); (c) 
\xi — X2\ < 0(l/n), \x[ — x'2\ < 0(1/71). These are just 
Yang's criteria 14] for superconductive ODLRO of p2 in 
the absence of ODLRO of pi. In case (c) p2 remains of 
order for arbitrarily large separation of the centers 
of mass X ^ {xi + X2)/2 and X' = {x[ + x'2)/2, the 
hallmark of ODLRO. On the other hand, in cases (a) 
and (b) p2 decays exponentially with |X — In the 
thermodynamic limit only configurations (c) contribute 
to the largest eigenvalue of p2, and p2 separates apart 
from negligible contributions (a) and (b) 13]: 

P2{xi,X2;x[,x'2) = n^sgn(a;i - a;2)e"^"l^i"^^l 
X sgn(a;'i - a;2)e~^"'^'i~^2l _|_ terms negligible for Ai. (3) 

By Yang's argument Q the largest eigenvalue is Ai = N, 
and this is confirmed by comparison with the Ai contri- 
bution Xiui{xi, X2)ui{x'i, x'2) to the spectral representa- 
tion of p2, implying that the corresponding eigenfunc- 
tion is ui{xi,X2) = C sgn(a;i — X2)e^'^"^^^~^^^ with 0] 
C = ^/n/L, confirming the value Ai — v? jC"^ — N. 
The range l/2n of ui is in the region of onset of a 
BEC-BCS crossover between tightly bound bosons and 
loosely bound Cooper pairs. There is an upper bound 
[1^ Ai < iV on the largest eigenvalue, so the untrapped 



FTG gas is maximally superconductive in the sense of 
Yang's ODLRO criterion. 

FTG gas on a ring: If the FTG gas is trapped on a 
circular loop of radius i?, with particle coordinates Xj 
measured around the circumference L = 27ri?, the FTG 
gas must satisfy periodic boundary conditions for both 
odd and even N because of single- valuedness of its wave 
function. Since the mapping function A{xi, ■ ■ ■ ,xn) = 
Y[i<j<e<N sgn(3^£ — Xj) is periodic (antiperiodic) for odd 
(even) TV as a result of its definition, it follows that the 
mapped ideal Bose gas used to solve the FTG problem 
must satisfy periodic (antiperiodic) boundary conditions 
for odd (even) N. The ground state of a FTG gas on 
a ring is then different depending on the particle num- 
ber parity. For odd N the FTG ground state in Eq. (P) 
is built from the zero-momentum orbital (j)o = 1/v^ 
and corresponds to mapping the FTG gas onto the ideal 
Bose gas ground state, the usual complete Bose-Einstein 
condensate (BEC), and is nondegenerate. On the other 
hand, for even N, which we henceforth assume, antiperi- 
odicity requires that the only plane-wave orbitals allowed 
are e''^^^ /vT with k = ^tt/L, ±37r/L, • • • . The ground 
state of this fictitious ideal Bose gas, and hence that 
of the mapped FTG gas, is then (N + l)-fold degen- 
erate, with energy eig envalue N{h^ /2m){Ti / Lf . These 
degenerate ground states are fragmented BECs with wN 
atoms in the orbital e^'^^'^^ and (1 — w)N in Q-^'^^i/L 
with < w < 1, and are conveniently labelled by a 
quantum number = [w — \)N = 0,±1,±2,- • • ,±y 
related to the eigenvalue P of circumferencial linear mo- 
mentum and that Lz of angular momentum z-component 
by P = £zh/R and Lz = Izfi. The angular momentum 
per particle is half-integral due to antiperiodicity of the 
orbitals, and the degenerate ground states are in one- 
one correspondence with the eigenstates of spin angular 
momentum z-component of N spin- 1/2 fermions. 

The ground state degeneracy makes the FTG gas on a 
ring a good candidate for detecting small external fields 
and linear accelerations. Suppose that there is a poten- 
tial gradient parallel to a diameter of the ring, or an 
acceleration leading to a gradient in the inertial poten- 
tial arising from Einstein's principle of equivalence, with 
the circumferential minimum of this potential occurring 
at a point xq. Then the degeneracy is lifted and to low- 
est order in degenerate perturbation theory all N atoms 
occupy the orbital (j^oix) = \/2/L cos[7r(a: — xo)/L], lead- 
ing to an observable asymmetric density profile n{x) = 
2n cos^[7r(x — xo)/L]. 

Due to its quantum coherence the FTG gas is also a 
good candidate for a sensitive rotation detector. Suppose 
that the ring trap is rotating with angular velocity uj per- 
pendicular to the plane of the ring. In the rotating coor- 
dinate system each atom sees an effective Coriolis force 
Fcor ~ 2mvxuj. Comparing this with the usual magnetic 
force -Fmag = {e/c)vxB, one sees that the kinetic energy 
operators in the Hamiltonian in the rotating system are 
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FIG. 1: Dependence of energies E on rotational flux $. Heavy 
line: Ground state energy Eo{^). Lighter lines: Lowest en- 
ergy for each value of total angular momentum. 



[pj - ^^]V2m where pj = {h/i)d/dxj, $ = ttR^uj is 
the Coriolis flux through the loop, and $o = h/2m is 
the Coriolis flux quantum. The energy of each state \£z) 
then becomes E = Eo{^ = 0) + M^iO^f ~ 24^] 
which is minimized when = if $ > and = — ^ 
if $ < 0, i.e., even a very small angular velocity leads 
to a nondegenerate ground state with all N atoms at 
either k ~ tt/L or k = — 7r/i. Generalizing to states 
differing from the $ = ground states by displacement 
in fc-space by integral multiples of I-k jL one obtains the 
^-dependent ground state energy £'o(<&) shown by the 
heavy line in Fig. ^ in which the lighter lines show the 
lowest energies for = ± -j, ■ ■ ' • The ground state 

energy is a periodic function of $ with period $0 in ac- 
cord with a general theorem fl3|, but unlike the usual 
situation for a superconductor, (a) there is no smaller 
period $o/2, and (b) for even iV, $0 = is a relative 
maximum of i?o rather than a minimum (as is the case 
of odd A^), the first minima occuring at $ = ±$o/2. 
The barrier heights of the energy landscape in Fig^ van- 
ish like for N 00, so flux quantization will not 
be observable for a macroscopic ring. However, it may 
be observable for mesoscopic rings using BEC-on-a-chip 
technology. For example, assuming a ring radius R =5 
/im, one finds that for ^Li, AE > ksT for T < 50 nK. 

Flow properties on a nonrotating ring: According to 
the FB mapping the excitation spectrum of the FTG gas 
is the same as that of an ideal Bose gas, and hence it is 
sufficient to analyze the latter. Since the excitation en- 
ergy of the ideal Bose gas is quadratic in the excitation 
momentum hq, the FTG gas does not satisfy the Landau- 
Bogoliubov criterion for superfluidity. We investigate 
here the possibility of flow metastability associated with 
barriers in the excitation energy landscape as a function 
of the transferred momentum. It was shown by F. Bloch 
[l6| that for the usual ideal Bose gas, which corresponds 
to the case of odd N in our treatment, no such barriers 



exist. In the case of even N both the ground state and the 
excitation branches are (iV -I- l)-fold degenerate, but it is 
sufficient here to consider the = ground state and the 
excitations arising from it by promoting atoms to higher 
fc-values. Generalizing Bloch's analysis, we note that for 
< < N/2 the lowest branch corresponds to excita- 
tion of V atoms from k = —'k/L to k = Zir/L, yielding 
a state with angular momentum z-component izti with 
tz — 2iy, and with excitation energy e{£z) — izfi^ /mR^. 
At I' = N/2 one has reached a state differing from the 
ground state by translation of all atoms by an amount 
27r/i in fc-space, and one can repeat this process, pro- 
moting atoms from fc — ir/L to Stt/L, yielding another 
straight-line segment connecting the points £z = N and 
£z = 2N on a parabofic curve {£zfi)^ /2NmR'^ , etc. To- 
gether with symmetry e{iz) = e(— ^z) this yields an ex- 
citation energy curve composed of straight-line segments 
as in the dashed curve of Bloch's Fig. 2 with the no- 
tation P = £zh/R. Hence for both odd and even N there 
are no energy barriers, and the FTG gas on a nonrotating 
ring does not exhibit flow metastability. 

Expansion from a longitudinal harmonic trap: We 
focus finally on a ID expansion, as could be achieved by 
keeping on the transverse confinement. If the ID interac- 
tions are suddenly turned off before the gas is let free to 
expand from a longitudinal harmonic trap, the density 
profile at long times reflects the initial momentum distri- 
bution . If instead the interactions are kept on during 
the expansion we find that the density profile expands 
self-similarly, while the momentum distribution evolves 
from an initial overall Lorentzian shape |12j to that of an 
ideal Bose gas. These properties can be demonstrated 
with the aid of an exact scaling transformation as 
we outline below. Since the FB mapping holds also 
for time-dependent phenomena induced by one-body 
external fields 17], the exact many-body wavefunc- 
tion tpp^xi. 
during the 
solution of 
tion for the 
an external 
with a;(0) = 
(j)o{x;t) = 
b{t) is the 



^N)]\^=l4>Q{.Xf,t) 



UJ^{t)b 
rt 



■ ■ ,XN;t) = A{xi, 

dynamics is fully determined by the 
the single-particle Schrodinger equa- 
orbital (j)Q{xj-,t). For the case of 
potential Vextix,t) — muj{t)'^x^ /2 
ujq the solution is known to be 

solution of the differential equation 
with 6(0) = 1 and 6(0) = 0, 



b 

T{t) = dt' 1/b^ and Eq ~ tiLLiQ/2. Since the unit anti- 
symmetric wavcfunction A is invariant under the scaling 
transformation, we immediately obtain the expres- 
sion for the many-body wavefunction, ^p(xi, ..,XN;t) = 

b-^/^i^Fixi/b, .., XN/b; 0)e^Wb^o^ Ef.i ^]/2xl.^-^NEorit)/H^ 

and for the one-body density matrix, pi{x,x';t) — 



; 0) exp 



This yields the mo- 
function of time. While 



mcntum distribution as a 
the intermediate-time dynamics has to be determined 
numerically, the stationary-phase method determines 
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FIG. 2: Momentum distributions of a FTG gas (solid lines) 
with A'^ = 9 particles as functions of the wavevector k at 
subsequent times t (in units of 1/ujo) during a ID expansion, 
and asymptotic long-time expression Q (dashed line). 

the long-time evolution of the momentum distribution 
in the same way as for the bosonic TG gas 0|. For 
the case of a ID expansion the scaling parameter is 
b{t) = ^y^+uJ^t^ and the momentum distribution tends 
to that of an ideal Bose gas under harmonic confinement, 

n{k,t ^ oo) c:i \ujQ/b\nB{kLUo/b), (4) 

where nsik) = 27rA^|^o(fc)Pi with ^o(fc) — 
'K'^/'^kosJ'^e''^^ /"^^'"'i^ and fcosc = l/a^osc- This be- 
havior is illustrated in Fig. |21 Quite noticeably, the 
"bosonization" time appears to be much longer than the 
"fermionization" time of the momentum distribution 
of the bosonic TG gas Note that the "dynamical 

bosonization" described above does not violate the 
Pauli exclusion principle: by using the above scaling 
solution for the one-body density matrix and fixing 
unit normalization of the natural orbitals at all times 
it follows that the eigenvalues aj of pi{x,x';t) are 
invariant during the expansion and hence always satisfy 
the condition aj < 1. 

In conclusion, we have found that (a) the un- 
trapped system exhibits superconductive ODLRO of the 
two-body density matrix p2 associated with its maxi- 
mal eigenvalue N and pair eigenfunction C sgn(a;i — 
X2)e~^"l^i~^^l ; (b) on a ring it has a highly degenerate 
ground state for even atom number, and it exhibits quan- 
tization of rotational Coriolis flux and high sensitivity to 
rotation and to accelerations, making it a good candi- 
date for high-sensitivity detectors; (c) the harmonically 
trapped system undergoes a "dynamical bosonization" of 
its momentum distribution during a ID expansion. 
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